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I Abstract 
< 

The critical and asymptotic behaviors of solutions of the sixth Painleve equation PVI , obtained in 
the framework of the monodromy preserving deformation method, and their explicit parametrization 
in terms of monodromy data, are synthetically reviewed and tabulated. 
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^ ■ 1 Introduction 

> ■ 

a ■ 

S In the last few decades, the Painleve equations have emerged as one of the central objects in pure mathe- 
matics and mathematical physics. They define new non-linear special functions, which find applications 
in a variety of problems, such as number theory, theory of analytic varieties (like Frobenius structures), 
random matrix theory, orthogonal polynomials, non linear evolutionary PDEs, combinatorial problems, 
etc. The properties of the classical special functions have been organized and tabulated in various clas- 
sical handbooks. The scientific community is now engaged in the project of a comparable organization 
and tabulation of the properties of the Painleve functions (this is part of the Painleve Project, see 
^T} " http://cio.nist.gov/esd/emaildir/lists/painleveproject/). 

In this paper, I provide two tables, one for the critical behaviors of solutions of the sixth Painleve 
equation PVI, and one for the formulas which parameterize, in terms of associated monodromy data, 
the integration constants determining the critical behaviors. 
• • ■ Though critical behaviors and expansions of the Painleve six functions at the critical points, and their 

parameterization in terms of monodromy data, have been extensively studied, the effort of collecting 
and tabulating them has been missing. It is the time to make a first point of the situation, in view of 
\ the Painleve Project. This is the aim of the paper. 

The tables are based on the results of Jimbo [17] and of my papers [10], [11], [12] and [9]. They 
represent essentially all have been obtained about the critical behaviors of the Painleve six functions 
and their parameterization in terms of associated monodromy data, by techniques of the monodromy 
preserving deformations method. 

The table of critical behaviors I propose here is essentially complete. Namely, the solutions which 
are tabulated form a set in correspondence to (essentially) the whole space of monodromy data of the 
isomonodromic Fuchsian system associated to PVI. 

Except for the case when logarithms appear in the expansions tabulated, the expansions are conver- 
gent and define true solutions. This follows from the results of [22] and [9]. When logarithms appear, no 
proof of convergence is known to the author, though one may expect that an extension of the method 
of [22] and [9] be possible for the power logarithmic expansions. 
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Bruno and Goryuchkina, in as series of papers which are summarized in their review [6], have com- 
puted the expansions of actual and possibly formal solutions of PVI. They have obtained all the expan- 
sions which can be constructed by means of their method of power geometry (see [4] and [5] ) . In section 
7, I compare their results with those obtained with the method of monodromy preserving deformations, 
showing that they coincide. 

The solutions (22) and (26) in the table of critical behaviors are also obtained by Shimomura [22], 
on the universal covering of a punctured neighborhood of the critical point. 



The sixth Painleve equation PVI is: 
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Following [7] (page 8), solving PVI means that: 

i) We determine the explicit critical behavior (or asymptotic expansion) of y(x), by an explicit 
formula in terms of two integration constants. 

ii) We solve the connection problem, namely we find the explicit relations among couples of integra- 
tion constants at different critical points. 

If we can solve i) and ii), then Painleve transcendents can be efficiently used in applications as 
special functions. The tables of the critical behaviors and their parameterization formulas in terms of 
monodromy data give answer to i) and ii). 2 

In the following, the subject "I" will be replaced by "we". Then, we recall the basic fact that, 
according to [18], PVI is the isomonodromy deformation equation of the 2x2 Fuchsian system below: 



~d~X 



= A(x,\)y, A(x,X):= 



A (x) , A x (x) Ax(x) 



A 



X-x A - 1 



As C. 



(i) 



where A(x, A) depends on x is such a way that 

A 12 (x,X) - 

Therefore, 

y{x) 



f(x){X-y(x)) 
A(A - 1)(A — x) 

x{A )i2 



x[(A ) 12 + (i4i)i 2 ] - (^1)12 (2) 

The 2x2 matrices Ai(x) depend on x in such a way that there exists a fundamental matrix solution 
$>(X,x) with monodromy independent of small deformations of x. They also depend algebraically on 
a, P, 7, 6 according to the following relations: 



A Q +A 1+ A x = -?fa 3 , 0^0. 



a=-(9 cc -l) 2 , 



Eigenvalues (Ai) 
1 



4 



0,1, x; 



7 = 



9 2 



2 1 2 x 



(3) 



(4) 



In the "A-plane" C\{0,x, 1} we fix a base point Ao and three loops, which are numbered in order 
1, 2, 3 according to a counter-clockwise order referred to Ao- We choose 0,x, 1 to be the order 1,2,3. 
We denote the loops by 70, j x , 71. See figure 1. The monodromy matrices of a fundamental solution 



2 A more restrictive definition of "solving" should include the determination of the distribution of the movable poles, 
for which we refer to [13] and [3]. 
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Figure 1: The ordered basis of loops 



*(A) w.r.t. this base of loops are denoted M , M x , M\. The loop at infinity will be 700 = 7o7x7i ; 
so Mqo = M\M x Mq. When A goes around a small loop around A = i, i = 0, x, 1, the fundamental 
solution transforms like i-> *Mj. Mo, M x , Mi are the monodromy matrices. There is a one-to- 
one correspondence between a point in the space of monodromy data and a branch y(x) (except when 
one Mi or M\M X M is the identity. Sec [11] and [12]). As a consequence, the integration constants 
of y(x) can be parameterized in terms of the monodromy data. Jimbo provided for the first time the 
critical behavior and its parameterization in terms of monodromy data for a wide class of two complex 
parameter solutions. His work [17] is the foundation of all subsequent developments. 

Organization of the paper 

Section 2. The table of critical behaviors is provided. 

Section 3. The table of the parametrization of basic critical behaviors in terms of monodromy data 
is provided. 

Section 4. The completeness of the table of the critical behaviors is discussed. 
Section 5. Okamoto's transformations are reviewed and the basic solutions arc defined. 
Section 6. The computation of critical behaviors through the matching procedure is reviewed. 
Section 7. The table of critical behaviors is discussed in details. A comparison with the results of [6] 
is provided. 

Section 8. Details on the parametrization of the critical behaviors in terms of monodromy data are 
given. The action of the Okamoto's transformations on the monodromy data is briefly reviewed. 
Section 9. The computation of monodromy data through the matching procedure is reviewed. 

2 Table of Critical Behaviors 

Here the table is given of the critical behaviors of branches of solutions of PVI at x = 0. The branch cut 
may be taken to be — 7r < arg x < ir. This is the first tabulation available in the literature. The detailed 
description and explanation of the critical behaviors tabulated below is in section 7. The numeration of 
the expansions in the table is that of section 7. The term "basic solutions" is explained in section 5. In 
order to understand the table, we introduce the set Hps — U ^351 where: 

T.% := [i y ^^+(-) k VT^2s)sgn{n(^2P + VT^25))}, k = l,2, (5) 
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or, in case \^2fi + {-) k \f\ - 26\ < 1, 

^ s :={[V^+(-) k VT^2s),-[^W+(-) k VT^25)}. (6) 

In the monodromy deformation parameters, this is Y,ps = {(0 O + 6* a; )sgn(3fJ(6'o + 9 x j), (9 — 9 x )sgn'3l((9o — 
6 X ))}, or = {(00 + 0*), (0 -0 X ), -(0o + 0x), -(0 -6 X )}. We also introduce the set n aj = Si^uO^, 
where 

^ 7 :={(V2^+(-) fe v ^)sgn(K(V2^+(-) fe v / 27))} J k=l,2 (7) 
or, in case |3£{\/2^ + {-) k ^2Fj}\ < 1, 

Q k ai := { (V^ + R fc v^r) , - (V^ + (-) fe V^) }• (8) 

Again, this is tt Q7 = {(0^ - 1 + ft) sgn(3f?(0 oo - 1 + ft)), (0c - 1 - 0i)sgn(3ft(0 co - 1 - 00)}, or 
- {(#co - 1 + 61), (600-I- -(0oo - 1 + ft), -(0oo - 1 - ft)} 
The behaviors in the table may be classified as follows: 

- Complex power behaviors (22), (26), (37) and (41). They are expanded in powers of x n+mX , for some 
n, m € Z and A € C. In this case, \y(x)\ may vanish, converge to a constant or diverge when x — > 0. 

- Inverse oscillatory behaviors (35) and (39). y(x) oscillates without vanishing when x — > 0, and may 
have poles in a sector centered at x = 0. 

- Integer power behaviors, namely Taylor series (30), (31), (45), (46), (27) and (38). 

- Power logarithmic behaviors, of type (28), (29), (33) and (43), namely power expansions with coefficients 
which are polinomyals of In a;. y(x) — > constant as x — > 0, where the constant may be zero or not. 

- Inverse power logarithmic behaviors of type (44) and (48), which can be expanded as asymptotic series 
of (Inx)- 1 . y(x) = 0(1/ lnx), or 0(1/ In 2 x), asi^O. 

Except for the case when logarithms appear in the expansions tabulated, the expansions are con- 
vergent and define true solutions. The proof is based on the results of [22] and [9]. When logarithms 
appear, no proof of convergence is known to the author, though one may expect that an extension of 
the method of [22] and [9] be possible for the power logarithmic expansions. 

In the table, a,<j> and a denote complex free parameters (integration constants), and v a real free 
parameter. The coefficients c nm , d nm and b n are rational functions of y/a, yf]3, ^ and \/T— 26. The 
coefficients b n (a) are rational functions of y/a, \//?, y/^y, \f\ — 26 and a. P n (lnx,a) are polynomials in 
In a; with coefficients which are rational functions of y/a, yf]3, ^7, a/1 — 2<5 and a. These coefficients can 
be recursively computed (essentially, by substitution into PVI. See [12]). In the case of (22), the c nm 's 
depend rationally also on a. In the case of (35), the c nm 's depend rationally also on v. 
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Complex power behaviors 


Free Parameters 


Other Conditions 


(22) y(x) = £~ i^E:=_„c„ ro (a^r 
Basic solutions 


IT, a^O 


0<9fcr<l, a&Xps, 
rn — 1 


(26) y(x) - £~ 1 E™=o ^K) m 
Basic solutions 




Rtr>-1, ff€Ej 7 , 
E^ 7 nZ = 0, fc = l,2. 

Cll = 1, 

c - - v ^ 


(37) y{x) = + E~i EJUo ^ ro (a^) m 


a 


Rp>-1, p + iefi^, 
Q* 7 nz = 0, fc = i,2. 
n£ 7 is (7) 


(41) - (1 + £~ 1 E„=o 4™(«^) m ) 


a 


a = 0. 
w g Z, sftcj > 0. 

we{v^,-V2^} if 7 < 0, 

otherwise w = sgn(5R v / 27)- 



Inverse Oscillatory Behaviors 


Free 
Param 


Other 


(35) y(x)= E~o^E™ + =-„-i^(e^ 2 -) ro ]" 1 
= [Asin(2i/lna; + 0) + J B + 0(a;)] _1 


u, 4> 


v e R\{o}, 2ii/ ^ o Q7 


(39) y(z) = [ 7ST g^ 7 = + ax^ + EZi * n E™=o c n+hm (ax^) m ] ' 


a 


z/eR\{0}, 2ii/e0^ 7 , 

fc = 1,2. 
ft* 7 c iR is (8) 



5 



Taylor expansions 


Free 
Par 


Other Conditions 


(27) y(*) = v 4^ x + ^ bnxn 




j/(ar) = if /3 = 0. 
V-2/3 ± VI - 2(5 g Z. 


(30) y(x) = El='i b n x n + ax^ + £~ w+2 M«)* n 
. . y(x) = a.x + a(a- 1) (7- a- i)x 2 + 

(31) +E^^(«K, " 7V ^° 


a 


One of V-2£ + VI - 2(5 or 

V-2/3 - Vi - 2(5 = TV e z 

and 

2/3 = -iV 2 , or 
\\f2a + V2a< - a/2T!- 
n {|7V|-l,l-|AT|}^0. 

(* 2/3 = 2(5 - 1 = when N = 0) 


(38) 2/W = ^ 5 +E"=i^" ; Basic Taylor 




V^i V27^Z 


(45) y(z) = E^o" 1 + + £~=|"l+i b n (a)x n 

(46) y(a;) = a + (1 - a)(<5 - + £~ =2 b n (a)x n , N = 

Two Basic Taylor solutions when N = 0, 1 


a 


One of V2a + V2~7 or 

V2Q 7 - V27 = n e z 

and 

2a = TV 2 , or 
{V-2/3 + VI - 2(5, V-2/3 - VI - 2(5} 

n {|jv|-i,i-|jv|} 9 6 

(* a = 7 = when iV = 0) 



G 



Power logarithmic behaviors 


Free 
Par. 


Other Conditions 


Ofl\ <ii(<r\ — vl-^l h <r- n 4- I n 4- h, , In <r 1 ^.1-^1+1 , 

(zsj y{x) — ij n=1 o n x + I a + 0|jv|+i ma; ix 1 + 
(29) y(x) = (a± ^-2/^lnx^x + Y^- 2 P n(^x;a)x n , ^ = 

Basic solution when N = 


a 


+ VI - 2(5 or 
y/-2P-y/l - 2(5 = N e Z 
(*2/3 = 25 - 1 if AT = 0) 

and 

V2^±^^ ±(\N\ - 1) 


(33) y(x)= [M±^ (a + i nx) 2 + _^] x + E - 2 p ri(lnx;a)x n 

Basic solution 


a 


2/3 ^ 2(5 - 1 


(43) j,(a;) = EL='o _1 ^ n + (a + b N lnx) X \ N \ +^ ={Nl+1 P n (lnx;a)x n 


a 


V2a + V27 or 

V2^-y/2^ = N e Z\{0}, 

and 

2a ^ TV 2 , 
V-2/3±\/l-2<5^±(|7V|-l) 



Inverse power logarithmic behaviors 


Free Parameters 


Other Conditions 


y(x) = {(a±V2^\nx)+J2n-i p n(^x;a)x n y 1 

(44) 

= ± 1 1 =F a +0 ( i 1 


a 


a = 7^0 


= ( [l_£ (a + l na; ) 2 + -^1 + J2n=i Pn+iilnx-a^Y 1 = 
(48) LL J J 

(7— a)ln' ! x L lna; Vln^~c/J 


a 


a ^ 7 



3 Table of Parametrization formulas 

The critical behaviors are governed by two integration constants, or free parameters, which appear in the 
explicit expansions at x = 0. We give a tabulation of the parametrization formulas of these constants, 
in terms of monodromy data of the associated system (1). As for their computation is concerned, the 
reader may see [17] and [2] for (22) and (26), and my papers [10], [11] and [12] for the other cases. 
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Paramctrization formulas are the achievement of the method of monodromy preserving deformations. 
They can be written at x = 0, 1, oo for the same monodromy data, namely for the same y(x). This allows 
to solve the connection problem. 

From the table of the critical behaviors, we extract a subset of behaviors which are basic (see section 
5 for a precise definition). All the other behaviors in the table can be derived from the basic ones, by 
the action of a birational Okamoto's transformation of section 5. All the behaviors at x = 1, oo can also 
be derived from those of the table, by an Okamoto's transformation (12) and (13), as it is explained in 
section 5. In this sense, all the information we need is contained in the basic behaviors listed below: 

o Solutions (22) and (26). 

o The power logarithmic behaviors (29) and (33): 

'(29) y(x) = (a±V^W^ x y + Y / n=2 P n(^x;a)x n , 2/3 = 25-1 

z + £~ =2 P„(lnz;a)z™, 2/3 ^ 25 - 1 

It is enough to compute the parametrization formulas of the integration constants for the basic critical 
behaviors (22), (26), (29) and (33). The other critical behaviors are obtained from the basic ones by 
birational transformations, and the parametrization formulas are correspondingly obtained by an action 
of the birational transformations on the monodromy data. Examples of these actions are formulas (52) 
and (53) to be introduced later, which allow to pass from x — to x — 1 and x = oo respectively, and 
(54), associated to the birational transformation (14). 

We give also the parametrization formulas for the solutions (35) and (39), which are not basic, but 
are important for their behavior and their applications in quantum cohomology (see for example [13]). 
We give also the parametrization formulas for the Taylor solutions below: 

(46) 1/ ( 3 ;)=a + (l-a)(^^ + ^ 2 ''r 1 («K, a = 7 = 

(47) y{x) = ±jfc + ax + Y%:Ma)x n , a ± 0, (v 7 ^ ± V^f? = 1, 2/3 = 25-1 

The latter of these is one of the two solutions (45) when N = 1. These solutions, together with (38) 
(which has no free parameter), are basic Taylor solutions, which generate all the other Taylor expansions 
via birational transformations. We will explain in section 5 that (46) and (47) can be obtained from 
(29). 

Let Pij =Tr(MiMj), i,j = 0,x, 1. 

Case of (22) and (26) . The two integration constants are a and a. Their parametrization is computed 
in [17] and [2]) (note that in [17], the last sign in formula (1.8) at the bottom of page 1141 is ±cr instead 
of to-)- 

1 fP0x\ 

a = —arcos , 

TT V 2 J 

= (Ox -do- <t)(0 x +0q- °-){0qo + 0i - a) 
° 4 ( x 2 (0 co +0 1 +a)F 

where 

= r(l + a) 2 r (1(0 O + 9 X - a) + 1) r (\(0 X - 0q - a) + l) > 
T{l-aYT{\{e + e x + a) + l)Y{\{6 x -e Q +a) + l) X 

y r (|(goo + 9i - a) + 1) r (\{9i - 0Qo - a) + 1) V 



(33) y(x) 



2/3+1-25 



(a + lnx) 2 + 



2/3 



2/3 + 1-26 
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and: 



U := x sin(7rtj)p a: i — cos(tt9 x ) cos^O^) — cos(ir9 ) cos(7T0i) 



+ — sin(7rcr)poi + cos(tt9 x ) cos(7r0i) + cos^floo) cos(7T0 o ) 

7T 7T 7T 7T 

V := 4 sin -(6» + X -<r)sm-(9 - 9 X + a) sm-(9 (x> + 9 1 - a) sin ^(#00 



01+ a). 



Note that for (26), one has to take the limit of the above a for a tending to one of ±(0o + @x) or 
±(#0 — Ox)- The limit exist finite and non zero. 



For (22) and (26), we have 
In particular 



Tr{M M x ) £ (-00, -2] U {2} 
Tr(M M x ) > 2 when cr e zR 



Case of solutions (35) and (39). The constants v and <f) of (35) in terms of the monodromy data are 
computed in my paper [12]. 

2cosh(27w) = -p 0x , ^ = * ln ( - ^) 
A is the coefficient multiplying the sine in (35): 



and 



where 



and: 



. 1 = ^ + (0oo - 1 - Kf¥* + O00-I + Ol) 2 } 



(0«> - 1 - ft - 2iv){6 tx> - 1 + 9 X + 2iv)(9 + 9 X + 1 + 2iv) 1 



8ii/(2iv -6 -9 x -l) 



, ._ r(i - 2iv) 2 r (1(0^ - 1 + ^ + 2w) + 1) r (fli + 1 - gpo + 2^) + 1) 
:_ r(i + 2ii/) 2 r (i(0oo - 1 + 0i - 2ii/) + 1) r + 1 - floo - 2^) + 1) 

r (|(go + 0. + 1 + 2^^) + 1) T (\{Q X - 80 + 1 + 2iiQ) V 

x r (±(0 O + flx + 1 - 2i V ) + 1) r - O + 1 - 2*1/)) w 
"1 



W := e 



27TJ/ 



• Sinh(27r^)pia; + 008(770^) COs(7r0oo) + COS(7T0 O ) COs(7T0i) 



+ 



1 



— - sinh(27T^)p i + 003(^0^) cos(7T0i) + cos(7T0 oo ) cos(7T0o) 



V := 4cos^(0 o + 9 X + 2w)cos^(0 o - 9 X - 2iv) sin^(0oo -1 + 6-1 + 2iv)sm^(6 00 - 1 - 61 - 2iv). 



As for the solution (39), we have 



2iv 



where r is the above. Note that in this case 2iv is equal to one of ±(0oo — 1 + 6\) or ±(0oo — 1 — 0i), 
and the limit of r above exists for these values of 2iv. 



For (35) and (39), we have 



Tr(M M x ) < -2 
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Case of (47): For the basic Taylor expansion (47), the parametrization is computed in my paper [10]. 
Note that 9^ ^ 1, one of Q$±6 X = 0, and one of 0^ — 1 ± 0\ = 1 or —1. The parameterization is: 

y{x) = 1 + a x + 0(x 2 ), 

1 — Poo 

9^ (2s + 9 x + l) 9 X [2 cosM^co + 6 X )) - p i] 

a = — : , s — 



2(floo-l) ' " 2[cos(7r(0 co -0 x ))-cos(7r(0 oo +0 x ))]' 
In this case 

Tr(M M x ) = 2, < M • M x , M Y > is reducible 

Case of (46): For the basic Taylor expansion (46), the parametrization is computed in my paper [10]. 
Note that one of 9^ — 1 ± 9\ = 0. The parameterization is: 

y{x) = a + X — 1(1 + 91-91) x + 0(x 2 ) 

1 Poi - 2cos(7r0 o ) 



1- S ' 4cos (£|il_^) cos (ZE| a + ^) 

This result can also be obtained from (9) below, via (14) and (54). In case of (46) we have: 

Tr(MoMr) = -2, < M • M x , Mi > is reducible 

Case of (29): The free parameter of (29) is computed in my paper [11], for 9 ca ±9i ^ 2Z. Note that 
one of 9 x ±9 a = 0. The parameterization is: 

a = e^ 9a 9 ^ E (9 a + l)+i7T + lE -^ - h " 



7T 
+ - 



4 sin tt0 o sin §(0^+01) sin §(0^-01) 
COS7t(6> + 0l) oj [cos7r(0 o + 0i ) - COS7t(0 O - @1 )] 



2 sin7T0 o sin |(0oo + Si) sin 1(6*00 — Si) 2 sin tt0 o sin 7r0i j 

Here is the Euler's psi-function, and 7# is the Euler's gamma constant. In this case 

Tr (M M X ) = 2, < M a ,M x > is reducible 

When 6 — > 0, namely /3 — > 0, (29) becomes the Taylor series (31) and 

2 COS 7T01 - poi 

° 4 S ill§(0 oo +0 1 )sin§(0 oo -0i) U 

Case (33): The free parameter for (33) is calculated in my paper [11], when 9q,9 x ,9i,9 qo £ Z. Note 
that 9 ± 9 X 7^ 0. The parameterization is: 

C - c + p 0x - p xl 29 x 
a= b=B + &t~&i +{lE ~ ) + 



2 2 / V2 2 y \2 2 / V 2 2 
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where 

4 

b = -(sinirfli s(0 o - x )s(9 + 9 X ) + sinTrflo 5(6*00 + 6> 1 )s(6» 00 - 6>i)) 

7T 



B 

2m 



c = 2 cos 7r(^o — 6*i ) 

1 



2 cos 7r(6>o + 6»i) + 4C0S7T61, costt^oo -4e i7rSl cos ir9 x - Ae-™ e ° costt^ + 3e" r ( 01 - e °) - e "( fl "- i) 

C* = 2e i7r<?1 cos7r6» x + 2 e - J7r8 ° costt^ - 2e" r(01 - e ° ) 
and s(z) := sin(7rz/2). In this case 

Tr(M M x ) = 2 

4 Completeness of the Table of Critical Behaviors 

The monodromy data 9q, 6\, 6 X and 6*00 are fixed (up to 6k i-> — #fc, fc = 0, x, 1, and 6*00 2 — 6*00) by 
PVI. Consider a triple of 2 x 2 matrices M , M x and Mi such that 

DetM = DetMi = DetM^ = 1, TrM^ = 2 cos 7^, = 0, 1, x, 00 

where :— M\M x Mq. These matrices arc defined up to conjugation by an invertible matrix. The set 
of such triples, modulo conjugation, is the set of possible monodromy matrices for the system (1) with 
the basis of loops ordered as figure 1. We call it, together with 9q, 9i, 9 X and 6*00, the monodromy data 
of Fuchsian systems (1), with the basis of loops ordered as figure 1. For an equivalent definition, see 
section 2 of [12]. 

To every branch y(x) of a solution of PVI, a system (1) is associated, and its monodromy group. 
Therefore, to y(x) a point in the space of monodromy data is associated. Conversely, to a point in the 
space of monodromy data, either one branch y{x) may be uniquely associated, or an entire family of 
solutions may be associated. Also, no solution of PVI may correspond, in principle, to some exceptional 
points in the space of monodromy data. In [11] we proved that a branch y(x) is uniquely associated to 
a point in the space of monodromy data if and only if: 

a*) none of the monodromy matrices M , M x , Mi and M\M X M (order of figure 1) is the identity. 

Let Pij :=Tr(MjMj), j = 0,x, 1 and p M :=TrM M = 2008 7^, fj, = 0,x, l,oo. It is known (see [16]) that 
Pox, Pxi, Pix and the p M 's (cquivalently, the M 's up to sign) are good coordinates for the monodromy 
data when: 

b*) the group generated by M , M x , Mi is irreducible. 

If a*) and b*) hold, we can parameterize the integration constants of the critical behavior of a branch 
y(x) in terms of the pij's and 6^'s, as in section 3. See section 8 for more details. 

A natural question is if the table of critical behaviors is complete, namely if it includes all the 
possible critical behaviors. A way to answer is to see if the correspondence {y(x)} {monodromy 
data} is surjective. Namely, the table is complete if the image, through this correspondence, is the whole 
space of monodromy data. 

The results of the previous section, supplemented by the application of birational transformations, 
bring to the following conclusion: 

To the set of complex power behaviors and inverse oscillatory behaviors (namely (22), (26), (37), 
(41), (35) and (39)) the whole subspace of monodromy data such that Tr(M M x ) ^ ±2 is associated. 
The subspace of monodromy data such that Ti(M M x ) = ±2 contains the image of all the remaining 
solutions in the table of critical behaviors (Taylor, power logarithmic and inverse power logarithmic). 
This image is the union of the following subsets (let N € Z): 

Tr(M M x ) = 2, 6 ± 9 X =/= and < M , M x ,Mi > contains no reducible subgroup 
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Tr(M M x ) = 2, one of 6 ± X = 2N and < M , M x > is reducible 
Tr(M M x ) = 2, one of 9 X - 1 ± 9 X = 2N + 1 and < M • M x , Mi > is reducible 

Tv(M M x ) = -2, - 1 ± 6i ^ and < M , M^Mi > contains no reducible subgroup 

Tr(M M x ) = -2, one of 6 ± 9 X = 2N + 1 and < M , M, > is reducible 
Tr(M M x ) = -2, one of 6^ - 1 ± 6i = 2N and < M ■ M x , Mi > is reducible 

In this sense, the table of the critcial behaviors is essentially complete. 

To conclude this section, we give more details about the above characterization of {y(x)} i-> {monodromy 
data}. 

- To the space of solutions (22) and (26), and their images through the birational transformation 
(14) (to be introduced in section 5), namely the solutions (37), (41), (35) and (39), there corresponds 
the subspace of monodromy data such that 

Tv(M M x ) + ±2 

- To the space of the power logarithmic solutions (28) and (29), there corresponds the subspace of 
monodromy data satisfying 

Tv(M M x ) = 2 cos 7T TV = ±2, one of 9 ±6 X = N, < M , M x > is reducible. 

Taylor solutions (30) and (31) are included in this case if moreover 

O = ±N, or one of 6^ - 1 ± 6>i = ±{\N\ - 1) 

- To the space of the power logarithmic solutions (33) there corresponds the subspace of monodromy 
data such that 

Ti(M M x ) = 2, 6 ±9 X ^0, < M , M x , M x > has no reducible subgroup. 

- To the space of the power logarithmic solutions (43), there corresponds the subspace of monodromy 
data satisfying 

Tr(M M x ) = -2 cosirN = ±2, N ^ 0, one of 6^ - 1 ± 6i = N, < M ■ M x , Mi > is reducible. 
The Taylor expansions (45) are contained in this case, when also 

600-1 = ±N, or one of 9 ±6 X = ±{\N\ - 1) 

- To the space of the Taylor solutions (46) there corresponds the subspace of monodromy data 
satisfying 

Tr(M M x ) = -2, 600 - 1 = Ox = 0, < M • M X ,M X > is reducible. 

- To the space of the inverse power logarithmic solutions (44) there corresponds the subspace of 
monodromy data such that 

Tr(M M x ) = -2, one of 600 - 1 ± 61 = 0, < M ■ M X ,M 1 > is reducible. 



12 



- To the space of the inverse power logarithmic solutions (48) there corresponds the subspace of the 
monodromy data such that 

Tr(M M x ) = -2, flco - 1 ± 0i ^ 0, < M , M x , Mi > has no reducible subgroup. 

- When 

< M ,M X ,M 1 > is reducible 

then, according to [14], all the corresponding solutions of (PVI) are equivalent by birational canonical 
transformations to the following one-parameter family of solutions, with 0^ + 0i + O + 9 X = 0: 

_ 0i + 0CQ - 1 + x(l + Ox) _ 1 x (1 - x) du(x;a) 
V[X) ~ 0^-1 0^-1 u(a:;a) dx ' Uj 

where u(x; a) = Ui(x) + au 2 (x); a e C, Mi(x) and ^(x) are linear independent solutions of the hyper- 
geometric equation: 

x(l ~ x)— + {[2 - (0^ + 0i)] - (4 - 9 X + 6 x )x} — - (2 - OO )(1 + 0> = 

Their behavior is again power logarithmic, Taylor or (37), but the parametrization in terms of po x , p x \, 
Poi does not apply. The rational solutions fall into the above case, see [20]. 

5 Okamoto's Birational Transformations and Basic Solutions 

Birational transformations are symmetries of (PVI), namely invertible transformations: 

y'^) = —/ x'^Mv (9o,9 x ,9 1 ,0 oo )^(9' o ,9' x ,0' 1 ,9' oo ) (U) 

such that y{x) satisfies (PVI) with coefficients , 9 X , 0i, 0^ and variable x, if and only if y'(x') satisfies 
(PVI) with coefficients 0' ,6' x , 9[,9' O0 an d variable x' . The functions P, Q are polynomials; p, <; are linear; 
the transformation of the M 's is an element of a linear representation of one of the following groups. 

* Permutation group: 

O = 1; e' x = 6 x , 0i = 0o, 0^=^co; y'(x') = l-y(x), x' = l-x. (12) 

9' X = 9 U 0i = x ; O = O , 0^ = 0^; y'(x') = ±y(x), x' = ±. (13) 

0o = ^-l, ^ = 0i, 0i = 0„ 0^ = 0o + l; y'(x') = ^, x = x'. (14) 

In terms of the coefficients of PVI, the transformation (14) is a' = —(3, (3' = —a, 7' = | — S, 5' = \ — 7. 

* Weyl Group of the root system D4: 

w\: 

0i = -0i; O = 9q, 9 x = X , 0^ = 0oo. 

a , 9 + 01 + 9 X + 0co , 00 + 01 — X — 0CX) 1 

y o - ^ ^ y i o 1" ^ 
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y x 2 ' 00 ~ 2 

0^ = 2-^; 0q = ^o, X = 6 X , 9' 1 — 6 1 . 

W4: 

O'oo = 2 — #co; 0' x = 2 — 9 X : 9' = 9q, 9[ =9\. 
The variable x' = a;, but y'(x) is quite complicated and will not be given here. 

* Shift lj :v= (vi,V2,v 3 ,v 4 ) ^ v + e,, j = 1,2,3,4, where e\ = (1,0,0,0), e 4 = (0,0,0, 1): 
h ■ 0' Q = 9 + 1, 9[ =^i + l; 6 X = 9 X , 9'^ = 9^. 

Q'o = Oq + 1, 9[ = 6\ — 1; 9 X = 9 X , 9'^ = 9^. 
9 X = 9 X + 1 , 9'^ = 6*co + 1 ; 9' = 9 , 6[ = 9\ . 
9 X = X + 1, 6>^ = 6*00 — 1 ; 9 = 9 , 6[ = 6\. 



h 
h 
k 



The variable x' — x, but y'(x) is quite complicated and will not be given here. 

Transformations w\, W2, W3, w 4 and I3 are sufficient to generate the other ij's, and they give a repre- 
sentation of the affine Weyl group of D4. The transformations of y(x) corresponding to w\,W2,Wz and 
W4 are obtained from formula (2.10) at page 355 of [21], while the transformation corresponding to Z3 is 
obtained from formulas (2.5), (2.6) at page 354 of [21], provided that h is substituted with h + of (1.13), 
page 352. 

The critical behaviors of y(x) when x — > 1 and x — >• 00 are obtained from the behaviors at x = 0. 
Actually, the transformation (12) exchanges the values and 1 of the independent variable. Therefore, 
when x — > then x' -4- 1. We obtain the behavior of y'(x') at x' = 1 from that of y(x) at x = 0. 
The transformation (13) exchanges the values x and 1 of the independent variable. Therefore, if x — > 0, 
.x' — > 00 and we obtain the behavior of j/'(x') from that of j/(x). 

The symmetry (14) allows to construct new behaviors at x = from known ones. An example are 
the inverse oscillatory behaviors or the inverse power logarithmic behaviors. See section 7. 

Definition: A set of basic solutions is a set of solutions that can generate all the other solutions of 
the equation via a birational transformation. In the same way, we call basic expansions those expansions 
that can generate, via a birational transformation, all the expansions which formally satisfy the equation. 

Birational transformations allow us to obtain all the solutions given in the table, form the set of basic 
solutions (22), (26), (29) and (33). 

In order to see this, first observe that the complex power behaviors, the Taylor expansions and the 
power logarithmic expansions given in the table are the only allowed forms of these kinds. This can be 
proved by direct substitution into PVI. 

As a second step, observe that the complex power behaviors (37) and (41), and the inverse oscillatory 
behaviors are obtained from (22) and (26) via (14). The inverse power logarithmic behaviors and the 
power logarithmic behaviors (43) are obtained from (28), (29) and (33) via (14). Details are explained 
in section 7. 

The third step is to prove that one can obtain all solutions (28) from (29). This is done by a birational 
transformation belonging to the affine Weil group of D 4 . This is because these trasformations can be 
used to change 9$ ± 9 X into 9o ± 9 X + N , N e Z. It is a formidable task to obtain y{x) by repeated 
applications of birational transformations, but one can avoid this by observing that a transformation 
(11) applied to a power logarithmic behavior necessarily produces a power logarithmic behavior or an 
inverse power logarithmic behavior. Then, the conclusion comes from the fact that the expansions in 
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the table exhaust all possible forms of such behaviors, being relations like 9^±9 V £ Z, or ^ Z, preserved 
by birational transformations of the afhnc D4. 

Finally, Taylor expansions can all be obtained from (29). For example, when we apply the birational 
transformation l 3 (or l^ 1 ) to (29), we obtain the case N = 1 of solution (28), namely: 

v 00 n— 1 

a - - ^2(3 (y=2j3 T l) (« - 7) In x J x 2 + ^ x n ^ b nm (a) ln m x, (15) 

' n=2 m=0 

^2/3 + VI - 2(5 = ±1 or ^2~j3 - Vl - 2(5 = ±1 
Then, let (3 = in (29) and a = 7 in (15). We obtain respectively the two basic Taylor solutions: 

y(x) =a+ ^—^-(2-/ -2a- l)x 2 + J2n=3 K{a)x n , this is (31) in the table 

y{x) = ±y/—2(3 x + ax 2 + b n (a)% n , this is (32), section 7 

These are transformed by (14) into (46) and (47), which may also be used as basic Taylor solutions. 
Then, birational transformations belonging to the affine Weil group of D 4 applied the basic Taylor give 
all the other Taylor solutions at x = 0. □ 



6 How to obtain the Critical Behaviors: Matching Procedure 

The matching procedure, supplemented by the symmetries of PVI, gives the critical behaviors at x — > 0, 
in the table. Originally, such procedure was suggested by Its and Novokshenov in [15], for the second and 
third Painleve equations. This approach was further developed and used by Kapaev, Kitaev, Andreev, 
and Vartanian (see for example [1]). It was implemented in my paper [10] for the case of PVI. We review 
here the method of [10]. 

Let us divide the A-plane into two domains. The "outside" domain is dehned for |A| > |x|' 5oc/T , 
Sout > 0. In this domain, (1) can be approximated by: 



dX 



A 



X 



xA x 



Nout 

E 

n=0 



+ 



A 1 



X 



^OUT, 



(16) 



The "inside" domain is defined for A comparable with x, namely |A| < |x|' 5/JY , Sjn > 0, and (1) can be 
approximated by: 



IN 



dX 



A a x 

X X- 



NlN 



71=0 



IN, 



(17) 



where Njn, Nout are suitable integers. 

The leading term of y{x) is obtained as a result of two facts: a) Systems (16) and (17) are isomon- 
odromic. b) Two fundamental matrix solutions ^/out(X, x), ^in(X,x) must match in the region of 
overlap, provided this is not empty: 



Vout(X,x) - *in(X,x), \x\ Sout < \X\ < \x\ SlN , x 







(18) 



This relation is to be intended in the sense that the leading terms of the local behavior of ^out and 
^ in for x — >• must be equal. It also requires that Sin < Sout- 



The simplest reduction is a Fuchsian systems, for Nj N = N ut = 0: 

d^OUT 



dX 



A + A x A x 



A 



A- 1 



OUT, 



(19) 



15 



IN 



dX 



Ao A x 
A A — x 



IN- 



(20) 



Practically, we match a fundamental solution of (16) for A — > 0, with a fundamental solution of the 
system obtained from (17) by the change of variables fj, := X/x, namely with a solution of: 



IN 



djjL 



^L + ^--xA 1 Vx>" 



IN, 



A 



(21) 



With the substitution A = £ in the system OUT, we can finally write a form for the matching which 
is well suitable in case that Nqut and/or Nin is geater than zero: 



dv 

d^iN 
d/j, 



Noui 



, ™ An + A, 



-xA x ]P (xuj 

n=0 

+ A ° + — 



1/(1/ -1) 



n=0 



The matching being 



fi fi(n - 1) 

*OE/t(^, x) ~ ^>in(h, x), v, fi -> oo 



IN 



br ijJV < |H < |a;|-* ot73 



|a;|-5oi7T-l < | M | < |a;|«/JV-l_ 



The above matching produces the leading terms of Ao(x), Ai(x) and A x (x) as x — > 0, therefore it 
produces the leading term(s) of y(x) in (2). 



7 Critical Behaviors when x — > 0. Description of the table 

We explain in detail the behaviors in the table. The behavior of case 1) below was first obtained by 
Jimbo in its original and fundamental paper [17]. The matching procedure, together with the Okamoto's 
transformations, allows to obtain all the other behaviors. 

In the paper [6] (and in the long series of papers of the authors of [6] which precedes [6]) the method 
of power geometry (based on Newton polygons) is used, in order to obtain all the formal expansions that 
satisfy PVI. They are classified in a way different from that we use here. In the course of the discussion 
below, we compare step by step the results obtained in the framework of the isomonodromy deformations 
method and the classification of [6], and show that they coincide. See the paragraphs denoted by an 
asterisk (*). 

1) Two-complex parameter (en, a e C) solutions [17], [9], [10], [11], [12], [2]: 



y(x) = E~ i % n YT m= - n c nm x m \ - 1 < ft<7 < 1, Cn^O 



(22) 



Note that the family is invariant for a ^ —a, therefore we can restrict to the range 



< 3?cr < 1 
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The c nm 's are rational functions of Cn and a, and algebraic functions of the coefficients of PVI. Note 
that in the table we have used a parameter a and put en = 1. Therefore, here, we have Cn = a of the 
table. If a does not belong to the set (6), then 



c i,-i = 



{V=2P- VT^2S) 2 -a 2 (7^2?+\/l — 25) 2 -a 



1 



16(7 4 



Cll 



(23) 



and c nm — q n m(V®i y/~$-> y/li Vl — 2(5; cr) x (cn) m , where g„ m are rational functions of their arguments. 
They are determined by the procedure of [12], section 7. If Sftcr > 0, the dominant term is 

y(x) = c 1 ,- 1 x 1 -°(l + 0(x°,x 1 -°)) 

On the other hand, if u = —2iv, i/eR, solution (22) becomes a three-real parameter solution 



n+l 



y(x) = x A sm(2u \nx + <j)) + B + ^ x n ^ c nm x n 



B 



B 



n—1 m=—n—l 

Av 2 + 2/3 + 1 - 2(5 



(24) 



Cnm = q nm (Va, y/jH, Vl, Vl - 2(5; v) x e im0 

where q nm are rational functions of their arguments, and the role of Cn as integration constant is played 
by (j) (which is just the log of cn). 

Remark: If a does not belong to the set (6) and 3ftcr > 0, we can also write (22) as 



y(x) = J2 x2n+1 E 



E- 2 " E 



ct, a„ 



G C 



n=0 



m— — n 



71=1 



m— — n 



oo 

Vr"(i + E : 



oo oo 



r m(l-p)+np\ ^ 



e c, w P ^ o 



(25) 



n—1 m—1 n=no 

where no = (m — l)/2 if m is odd, no = m/2 — 1 if m is even. 

The leading term of the expansions (22) was obtained first by Jimbo in [17]. It can also be obtained 
with the matching procedure of [10], making use of IN and OUT Fuchsian reductions of the system (1). 
The full expansion (22) can also be obtained by local analysis, with the method of [22], [9]. This method 
proves also the convergence of the expansions, which therefore define true solutions of PVI, which are 
basic solutions. In [12] a recursive procedure is given to compute the coefficients c nm . 

(*) The class of the solutions (22) when $ta ^ is the class Ao of [6]. In [6] the _4 -solutions are 
written as (25), though the range of n, m is looser, namely n + m > 0, n, rn > 0. 

(*) The image of (B^ U Bq) U (B^f U B^ ) of [6] through the symmetry (14) is in the sub class of 
solutions (24). 



2) Suppose that PVI is such that y^2(3 + Vl - 2(5 £ Z or yf^Zfi - Vl - 2(5 £ Z (namely, 6 + 0, 



or 9 — 9 X £ Z). In this case, PVI admits a one-complex parameter (cn G C) family of solutions (see 
[10] and [12]) 



y(x) = E~ i x n E^ =0 c nm x m ° = Er=o Vk(x) (c n x°)\ a G S^, 9for > -1 



(26) 
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The expansion is convergent, as it follows from the local analysis of [22] and [9], so it defines true solutions 
of PVI, which are basic. The coefficient Cn is the free parameter, en = a of the table. The c nm are 

functions q nm (\/ci, VP, y/l, VT~~2S) x (cn) m , like in the case of (22). In particular cio = 



-2/3 



There are two possible solutions, corresponding to the signs ± above, namely to the fact that a € g , 
k = 1,2. We stressed that a is not a free integration constant. The functions yu{x) are convergent 
Taylor expansions of the form 

y Q (x) = ^ / _ + V ( _^^___ ? a; + J2 gno(V^, V^, V7, 7T^5)x", y (x)=0if/3 = (27) 

j/, (a:) = O(^). 
If — 1 < jRcr < 1, the leading terms of y(x) are 

y(x) = ^ ^ x + cnx 1+a + ... 

If (5 = 0, the leading term is always cui 1+,T and 

y(x) = Cll x 1+ff + ]T x" ^ c 



oo n 

mo- 



n— 2 m—1 



The solution (26) reduces to a Taylor series solution without free parameters, when en = 0, namely 

y(x) = y {x) 

The method of local analysis of [22] and [9] allows to extend the result of [17] in the case jftcr > 1, 
a € and this precisely yields existence and convergence of solutions (26). Note that when —1 < 
jfter < 1, they are a subclass of (22). Actually, the coefficients of (22) vanish when a belongs to 

the set (6) (see formula (23)) provided that \/—2(3 ± y/1 — 25 Z (otherwise coefficients diverge). Also 
note that this happens when we choose Cn ^ as integration constant. One can choose c\-\ instead, 
and express Cn from (23). In this case, again the coefficients of xl m l CT vanish for a in the set (6) and we 
again obtain (26) with a ^ —<r, ^St(-a) > — 1. 

(*) The symmetry (14) transforms the set (Bi(non log) U B 2 (non log)) U (B^ U Bf) U (B% U B%) U 
B 6 U (Bg~ U Bg~) of [6] into the solutions (26) with /3 ^ 0. It transforms Cg° U B 8 U B 9 of [6] into the 
solutions (26) with /? = 0. This will be clear after reading 2.1) below. [Here we use B^(non log) to 
denote the subclass of Bi such that the coefficients of powers of x are constants, and the power of x have 
no integer exponent (see 2.1) below). We distinguish it from the subclass Bj(log) whose coefficients are 
polynomials of In x, and the powers of x have integer exponents (for some value of the coefficients of PVI 
the coefficients may be constants, and the solutions become Taylor series. See 3.1) below). Bi = B^(non 
log)U^(log).] 

We conclude that the class A U (B$ U Bq) U (Bi (non log) U B 2 (non log)) U (S; 1 " U Bf ) U (B% US 2 ")u 
B 6 U (Bq U Bq) U (Bf U Bf ) U B 8 U B 9 U C£° of [6] is contained in the union of 1) U 2) U {the image of 
1) U 2) through the symmetry (14)}. 



3) If + VI - 2(5 = N or - VT" 7 25 = N (namely, O + X = N or 6 - 6 X = N), 



N G Z, then solution (26) is no longer defined, because some coefficients diverge. We have instead 
solutions formally like (26), where a is the integer N and the divergent coefficients c nm are replaced by 
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polynomials of lnrr. Namely, from the results of [10], [11], we see that there exist one-complex parameter 
(a G C) solutions 



y(x) = Y}^b n x n + (a + b m+1 lnx)^ 1 + S~ |JV|+2 P n (lna;; a)x n , N ^ 



(28) 



y(x) = (a±y/-2l3\nxJx + J2n=2 P n(lnx-,a)x n , 2 ^ N = 2 ^°_ 1 



(29) 



where 6„ are rational functions of y/a, \f]3, y/j, and y/1 — 26; a is the free complex parameter; P„(lnx; a) 



are polynomials of In x of degree n— \N\ , with coefficients that are rational functions of \fa, \f]3, ^7, \/l — 28 
and a. Solution (29) is a basic solution. In particular, the above (28) and (29) reduce to convergent 



Taylor series (studied in [10] and [19]) if either {\/2c^ + ^2% V2a - n {\N\ - 1, 1 - |iV|> ^ 



namely + 6>i - 1, 6^ - 0i - 1} n {|iV| - 1, 1 - \N\} ^ 0, or if 2/3 = -N 2 , namely O G {-/V, -iV}. 
The series, starting with ir-term, are 



y(x) = El2i &nX n + ax^ 1 + E~ ijvi+2 &n(a)x", N ^ 



(30) 



= y(x) = ax + ^(2 7 - 2a - l)x 2 + E^ =3 b n {a)x n , ^ = ^ _° = Q 



(31) 

Here 6„(a) are rational functions of \/a, \/]3 7 ^7, \/l — 2<5 and a. In (31) we have /3 = 1 — 25 = 0, namely 
9 x = e = 0. It is the limit of (29) for /3 ->• 0. 

In [10], the solution (29) is obtained with the matching procedure, making use of Fuchsian IN and 
OUT systems as reductions of (1). Convergence of (28) and (29) is not proved, but we expect it, the 
proof being possibly based on an implementation of the method of [22] and [9]. All the coefficients can 
be recursively determined by direct substitution into PVI. 

The Taylor series are obtained by non-Fuchsian IN and OUT systems as reductions of (1). In [10], 
three basic Taylor series are obtained in this way, with non vanishing constant term . They will be given 
later, see the solutions (38), (46) and (47). The symmetry (14) transforms (46) and (47) respectively 
into (31) and the following: 

00 

y(x) = x/^2f3x +a x 2 + ^b n {a; ^hl, yfp)x n , (3^0, ± VT^ 2S) 2 = 1, a = 7 (32) 

n— 3 

The last solution is included in (30) for N = 1. Taylor expansions are also computed in [19], and their 
convergence is there motivated. 

Solutions (31), (32) and (27) are basic Taylor solutions, namely all the other Taylor expanded so- 
lutions at x = can be obtained from them by birational transformation. For example, when N = 1, 
solutions (30) contain two nonequivalent families, namely the family (32) and the following family of 
solutions: ^ ^ 

y(x) = x + ax + -<z(7 - a)x 3 + 0(x 4 ), S — -/3 = - 
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The latter is equivalent to (31) via birational transformation, while (32) is not. 

(*) The symmetry (14) transforms the class Si (log) U B 2 (log) of [6] into (28), and B A U B 5 into (29). 



4) If 2/3 7^ 25 — 1 (namely, 9o±0 x ^ 0), besides log-solutions (28), there are one-complex parameter 



(a e C) logarithmic solutions ( see [17], [10] and [11]), which are also a basic solutions: 



y(x) 



2/3+1-25 
4 



a + lnx) 2 + 



2/3 



2/3+1-25 



X + J2n>2 p n0nx;a)x n , aeC. 



(33) 



The above solutions are obtained in [10] or [11] by the matching procedure, making use of IN and OUT 
Fuchsian reductions of system (1). As far as convergence is concerned, no proof is published, though we 
expect that it is convergent and an implementation of the method of [22] should allow to prove this. All 
the coefficients can be recursively determined by direct substitution into PVI. 

(*) The symmetry (14) transforms the class B3 of [6] into (33). 



I ■ the following, we apply the symmetry (14) to the solutions l)-4), and we obtain new solutions. In 
all which follows, it is understood that wherever a, (3, 7 and S appear in formulas related to l)-4) (for 
example, the coefficients of the expansions), we must substitute (a, /3, 7, 6) M> (— /3, —a, | — 5, | — 7). 



1.1a) When a does not belong to the set (6) and > 0, then from (22) we obtain the following 
solutions, via the symmetry (14): 



y{x) 



1 



x a /ci, 



00 00 

Ci _i ^— ' ^-^ 

' n— m——n 



n+2 e n + i, m _i mg 
" ci,_i X 



Namely 



y(x) = ^ a;" ^ 4 m i mp , p = 1 - ct, < < 1, d nm e C (34) 

n— 1 m——n 

Note that 1 + p does not belong to the set (8). If we do the explicit computation above, we find that 

d nm = for in > n 

Thus, (34) is a solution (22). 

(*) We have re-obtained the class Ao of [6]. 

Note: To verify that d nm = for m > n, one can do the explicit computation of the d nm by ex- 
panding x/y(x) and computing the coefficients. Note that to get the correct terms in the truncation 
Y,n=i xn Y,m=- n dnmX mp , one needs to start from (22) truncated at J2n=i a; ™ EL-» c nm x mp . One 
can also conclude directly from the form (34) that d nm — for m > n. Actually, such a form, substituted 
into (PVI), implies d nm = for m > n. Namely, we always have that: 



y( x ) = Yl x " rf " 



rap 



substitution into (PVI) 



y 



n—1 m— — oo 



n—1 m——n 
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1.1b) From (24) via the symmetry (14), we obtain a three-real parameters (v £ R, (f> £ C) solution 
(see [12]): 



y(x) 



i -1 



1 



A sin(2y In x+^)+B+0(a:) 



(35) 



where 



rs+l 



im<f> 2imu 
X , 



n— 1 m— — n— 1 



The expansion is convergent. All c nm , are rational functions of u, and ^/a, \//?, y/j, vl — 2<5. If argx 
is bounded (namely, when we considered y(x) as a branch), then (35) admits two infinite sequences of 
movable poles in the neighborhood of x — 0, asymptotically distributed along two rays and accumulating 
at x = (see [13]). One can do a formal (non convergent) expansion, when x is far from the poles, as 



1 



oo oo 



^sin(2z^lnx + <j>) + B 



+ 5 v x 2iv \ l + ^>" S « 



8v, Snm € C (36) 



n— m=— n 



(*) The class of solutions (35) coincide with the union (B£ UBq) U (B^UBf) of [6], where a formal 
expansion, in a form equivalent to (36), is given. 



2.1) The symmetry (14) applied to case 2) gives one parameter (a £ C) solutions, defined when at 
least one of 



2a + ^27 £ Z 



We distinguish two cases, a ^ and a = 0. 



2a - V27 g Z (namely 6>oo +0i ^ Z or 0oo - #i £ Z). Let w G f2 a7 - 



o Case 



2.1a) When $ta > in (26) we obtain the following, via the symmetry transformation (14): 

1 



y{x) = 



cio (1 + + ET=i £™= Z^x* 



OO OO 



n— ?7i—0 



where 601 is the arbitrary constants, and the other coefficients are rational functions of yfa, \f$, ^7, \/l — 25 
and 601 • In particular 600 = • With the substitution uj = p + 1, the above becomes 



(37) 



p + l£^ 7 given by (7), fc = 1, 2, => 3ft/? > -1 
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Here a is the arbitrary constant, d nm e C are rational functions of y/a, \fft, y/j, and \f\ — 2(5. The 
yi(x) = 0(x l ys are Taylor series. In particular 



(38) 



Note that yi(x) can be directly obtained from (27) by the symmetry (14). Its convergence is proved in 
[19]. 

(*) In terms of the classification of [6], the class of solutions (37) coincides with the union of the 
class Bq with the subclass of B\ U £>2 which does not contain logarithmic coefficients. 



2.1b) When 3?cr = in (26), we obtain through (14) 



(39) 



Here 2iv = uj £ Q a 0- Note that Q a p must contain an imaginary element. We have two solutions above 
with lo = ± (v / 2a~+ v / 27) , when \f2a+yf2rf is imaginary, or we have two solutions with u> = ±{\/2a— \f2q) , 
when \/2a — y/2rj is imaginary. In particular 



to(x) 



+ ax 



liv 



— = — ^t—Ti — = + {cos(2^1na; + d>) + isin(2^1na; + d>)\, 
V^+(-) fe V7 1 J 



<t> = In a, 



1iv£ Q k aj = {(V2^+(-) fe V2^),-(\/2^+(-) fe V2^)} CiR, fc= 1,2. 



The coefficients c n + i. m 's are rational functions of i/a, ypy and \l\ — 2b. The denominator X)^°=o a; ™ T n( :z ') 
may vanish, so y(x) may have movable poles. In a neighborhood of ir = they are asymptotically dis- 
tributed along the ray of the zeros of tq(x). The latter are given by 



\x\ = ex P j^; ( ar S~ + ( 2l+ 1 ) 7r )} ' arg x = -^; ln 



cio 



ClO 



Far from the line of the poles, we can write 

1 



, /€Z, i<0 



y{x) 



tq(x) L r (a;) 



n(x) 
To(ar) 



= 0(z) 



(40) 



(*) The class of solution (39), when formally expanded from (40), coincides with the union (B+ U 
Br) U (B+ U B 2 ") U (B+ U fig" ) of [6]. 



-2.1c) When —1 < 3?ct < in (26), we obtain a solution for —1 < JJcj < of the form 

y{x) = 



1 | cio oj i y^oo „ v^n c„ + i, m 
1T cii ^Z^n=l x Z^ TO =-1 en 
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CO CO 



CO CO 



ax 



^ ^ x ^ ^ b nrn x — ^ x ^ ^ d nra x ^ 



n—0 m——n 



n— 1 m— — 7i 



where a = 1/cn is an arbitrary constant. This should coincide with a solution (22) with a = 1 + cj, 
< 9ftcr < 1, namely 



ax 



n—l m——n 

(*) In terms of the classification of [6], we have obtained solutions in Aq- 



o 

o Case 

o 



a = 



2. Id) The symmetry (14) applied to case 2) when (3 = yields the case a = 0: 



y(x) = 



Namely 



(41) 



w - V^7 sgn(3?V27) or u e { - V^tI if -l<»V^f< 1 ; 

Here a is the free parameter and > —1. d nm are rational functions of y/a, \/]3, ^ and \/l — 25. 
yi{x) = 0(x l ) are Taylor series, and yo(x) is monic, namely yo(x) = l + 0(x). The above \y(x)\ diverges 
for x — > when 3?cj > 0. When g iR, it becomes the special case a = of (39). 

Note that when — 1 < SRu; < 0, we obtain solutions of type (22). On the other hand, we obtain truly 
new behaviors (divergent, or oscillatory and non vanishing) when 



oj = ^27 sgn(3?y2x) g Z, or we {^2^, -^2^} C iR when 7 < 



(42) 



(*) The class of solution (41), with restriction (42), coincides with C§° U Bg U B$ of [6]. 



3.1) The symmetry (14) applied to case 3), gives solutions existing when \[2a. + y/2j = N 



or 



y/2a - V27 = N (namely 6^ - 1 + 6 X = N, or 6 X - 1 - 6>i = TV), AT e Z: 



(43) 
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(44) 



The 6„'s are rational functions of y/a, V/3, V7 and y/l — 26; a is free complex parameter; P„(lnx; a) are 
polynomials of In x of degree n— | JVj+1, with coefficients which are rational functions of ^fa, V/3, VX \fT~~25 
and a. The logarithmic terms vanish, and the solutions become convergent Taylor expansions, either 

if 



(45) 



W-2P + VI - 26, V-20 - VI - 26} n {\N\ - 1, 1 - |iV|} ^ 


, namely {6» +0 X ,0 O - #*} n {|JV| 


1 - |JV|} ^ 0, or if 


2a = N 2 


(namely 9^ = N + 1). The expansions are 



= Lilt 1 ^ n + ax^l + ££=1*1+1 b n (a)x n , N ^ 



y {x) = a+(l-a)(5- f3)x + En=2 b n(a)x n , N = 0, a = 7 = 



(46) 



We observe that the case (45) with N = 1 gives two subcases: 



j/(a;) = — L + ax + 0(x 2 ), a ^ 0, (V2a± V^t) 2 = ^ 1-25 + 2/3 = (47) 
V2a 



and 



y(x) = l + ax+ia(l+4a + 2/3-2<5)x 2 + 0(x 3 ), a = ^, 7 = 



The latter is equivalent to (46) via a birational transformation. 

Solutions (46), (47) and (38) are basic Taylor solutions, which generate all the other Taylor expansions 
by the action of the birational transformations. They are obtained in [10] making use of the matching 
procedure with non Fuchsian IN and OUT reductions of (1). They are also derived in [19]. 

(*) In terms of the classification of [6], the class of the solutions (43) is the subclass of B\ U Bi with 
logarithmic coefficients. The class of the solutions (44) coincides with the class B4 U B 5 . The sub case 
(46) gives the class B w of [6]. 



4.1) Symmetry (14) applied to 4) gives the following solutions defined for a ^ 7 (namely $^±6*1 7^ 



1): 



(*) In terms of the classification of [6], this is the class Bz- 



(48) 



(*) Final comparison with the results of paper [6]. We have shown that the set of the behaviors 
at x = obtained in the framework of the isomondormy deformation method, in particular using the 
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method of [17] and the matching procedure of [10], coincides with the set of expansions at x — obtained 
by power geometry in [6]. Summarizing, at x = 0, we have: 

Behaviors of [17], [22], [9], [10], [11], [12] 39 families of [6] 

1) behavior (22), with 3£er ^ 
[ This includes cases 1.1a), = A® 

2) with -1 < sftcr < 1, and 2.1c)] 

1.1b) behavior (35), and 1) behavior (24) = U B$ U U £f and image 1 

2.1a) behavior (37), and 2) behavior (26) if Sftcr > = ,8i(non log) U ^(non log) U Bq and image 

2.1b) behavior (39), and 2) behavior (26) if Via = = B\ U £f U 6 2 + U U B 6 + U and image 

2. Id) behavior (41), and 2) behavior (26) if /3 = = Cg° U B$ U £ 9 and image 

3.1), and 3) = Si (log) U B 2 (log) UB 4 U B 5 U Si and image 

4.1), and 4) = B3 and image 

1 "image" means the image of the families indicated through (14). 

In [6] a classification of actual, asymptotic or formal expansions which satisfy an ODE is given. It terms 
of it, we have that: 

- Solutions (22) and (26) for 5ftcr ^ 0, (37) for 3fy ^ 0, (41) for ^ 0, and the Taylor series (30), (31), 
(45), (46), (27) and (38) have power expansions. 

- Solutions (28), (29), (33) and (43) have power logarithmic expansions. 

- Solutions (44) and (48), when expanded in powers of (In a;) -1 , have complicated expansions. 

- Solutions (24) (namely, (22) for Kcr = 0) and (26) for sftcr = 0, (37) for 3£p = 0, and (41) for Ku; = 
have semi exotic expansions. 

- Solutions (35) and (39), when formally expanded (like in (36)), have exotic expansions. 



8 Details on the Parameterization in terms of monodromy data. 
Connection Problem 

As already esplained, to every branch y(x) of a solution of PVI, a system (1) is associated, and thus 
to y(x) the monodromy group of the system is associated, namely, monodromy data are associated 
(definition of monodromy data is in section 4) . 

It may well happen that to a family of solutions y(x) the same monodromy group is associated. An 
example, which is found in [11], is a one parameter a e C class of solutions existing when 

a = 7 = 0, 5 = i /3 = -2p 2 ; P eZ, p^0 

Their critical behavior is power logarithmic or inverse power logarithmic at the three critical points: 

x [l -p 2 (\nx + a) 2 ] , x -> 0, 

y(x)~t 1 - p 2 (In i + poo) 2 , x ^00, (49) 



1 



7: 



1. 
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where: 



7r(41n2 - 1 + a) 



7T 



.2 



Poo = 



7r-i(41n2- 1 + a) 



21n2 + l, pi = 



41n2 -1 + a 



-ln2 + l. 



The monodromy group associated to (49) is independent of the parameter a, and is generated by: 



In particular, Tt(M M x ) = Tr(M Mi) = 2 and Tr(MiM x ) = -2. 

Cases like the above example occur for special values of the coefficients of PVI. In general to y(x) 
is associated a point in the space of monodromy data and to this point no other branch than y(x) 
corresponds. In section 4 we have seen that this occurs if and only if none of M , M x , Mi and Mx, is 
the identity. Moreover the values po, p x , pi, p 00l p$ Xl p x \ and poi are good coordinates of the subspace 
of monodromy data such that the group < M , M x , Mi > is not reducible and none of the matrices M , 
Mi, M x and Mx) is the identity. 

Let the branch y(x) depend on (at most) two integration constants c\ and c 2 : 



for example, in case of (22), (ci,c 2 ) = (a, en) (or (a, a) in the table). When none of Mo, M x , Mi and 
Moo is the identity, and < Mo, M x , Mi > is not reducible, y(x) is uniquelly associated to a point in the 
space of monodromy data, and the integration constants ci,ci are uniquely determined functions: 



[Note: (p ,p x ,pi,p oo ), or equivalently (0 O , ), are given by the coefficients of PVI. Only two of 

POx,Pxi,Poi arc independent, because they are linked by the relation MiM x M = Moo, which is a cubic 
in terms of the p^ and Pij. Thus, two are the genuine free complex parameters]. 

In case Ci,C2 appear as the constants in the explicit formula of the critical behavior of y(x) at one 
critical point, then the method of monodromy preserving deformations allows to compute (50) explicitely. 
For a given point in the space of monodromy data, namely for a given y(x), the parameterization can 
be done at each of the critical points x = 0, 1, oo. In this way the connection problem is solved. 

In order to compute the monodromy data of (1) associated to a critical behavior of y(x) at x = 0, 
we procceed as follows. Once the matching ^out ^ ^in in (18) has been completed - and thus the 
critical behavior of y(x) has been determined - we have to match ^out with a fundamental solution \E" 
of (1) at A = oo, and we have to match ^in with the same '5 in another region of the A-plane, typically 
around A = or x. If this is done, then Mi of If coincides with Mi UT of ^of/T, while Mo and M x 
coincide with M$ N and M* N of */jv- 

The crucial point is that we are able to compute the monodromy matrices M° UT , Mq N and M^ N 
of the reduced OUT and IN-sytems exactly, namely that we are able to solve these systems in terms of 
linear special functions. 

The entries of the monodromy matrices so computed, depend on the integration constants c\ , C2, 
namely they depend on the free parameters which govern the critical behavior of y(x). In case y(x) is 
uniquely associated to the monodromy matrices just computed (namely, none of Mo, M x , Mi and Moo 
is the identity, and < M , M x , Mi > is not reducible), we calculate the traces of MiMj 



[ P01 — P0l(ci, C2, Oo, X , ffi, Ooo) 

These are classical functions of (a, C2, 0o,0 x ,0i,6oo). This follows from the fact that they are obtained 
by means of asymptotic reduction of (1) to the systems (17) and (16), which can be solved in terms of 




y(x) = y(x,a,c 2 ) 



Ci(P0,Px,Pl,Poo,P0x,Pxl,P0l), i = 1,2 



(50) 




(51) 
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classical functions (for definition of classical functions see [23], [24] and [25], or the introduction of [12]) 
The system (51) can be solved in c\ and C2, and we obtain (50). Moreover, the functions in (50) are 
classical of their arguments. They are the parametrization formulas of section 3. 

It is enough to give the parameterization at x — 0, because the critical behaviors and the parame- 
terizations at x = 1 and oo can be obtained via (12) and (13) provided we change the parameters in the 
parameterization formulas according to 



'Poi 



POx 



Pxl 



for (12), and 

' POx = 

Poi = 



P'u 



-POI - POxPxl + PooPx + PlPQ 
Pxl Ol < 

Pax 

-POI - POxPxl + PooPx + POPl 

Pox or '- 

Plx 



' Poi = 


-Poi 


POxPxl 


+ P'ooPx + PlP'o 


< POx = 


P'lx 






- Pxl = 


POx 






' Pm = 


—POx 


- P'oiPxl 


+ P'ooPl + P'aP'x 


< POx = 


Poi 






> Plx = 


P'lx 







(52) 



(53) 



for (13). 

Note also that it is enough to compute the parametrization formulas for the cases N), N=l,2,3,4, of 
section 7. One can obtain the parameterization for the solutions in cases N.l), from that of the equivalent 
cases N) [equivalent via symmetry (14)] . As it is explained in [12], section 4, one needs to make the 
substitution 

(jPo,Px,Pi,Poo;Pox,Pm,Pxi) i-> (-Poo,Pi,Px,-Po;-Pox,-Poi,Pxi), a ^ 1 - a (54) 
in the formulas of cases N) . 



9 Further Details on the Matching to obtain Monodromy Data 

The standard choice of a fundamental solution of (1) is 

[/ + 0(±)] A-^ ff3 A fl ~, A^oo; 

Mx) [I + 0(A)] \^\ R °C , A -> 0; 

^(a0[J + O(A-a0] {\-x)^{\-x) R °C x , A -> x; 



*(A) = < 



(55) 



k Mx)[l + 0(X-1)] (A -1)^-3( A _ \)RiC u 



A4l; 



Here tpo(x), ip x {x), ipi(x) are the diagonalizing matrices of Aq(x), A\(x), A x (x) respectively. They are 
defined by multiplication to the right by arbitrary diagonal matrices, possibly depending on x. C v , 
v = oo, 0, x, 1, are invertiblc connection matrices, independent of x [18]. Each R v , v = oo, 0, x, 1, is also 
independent of x, and: 



i?„ = if 9 V <£ Z, i?^ 







f(o 


I) 








o) 



if 0^ > integer 
if ^ < integer 
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If 6i = 0, i = 0,x, 1, then Ri is to be considered the Jordan form ^jj ^ J of Ai. 

(*) Matching * ^ $>out- 

A = oo is a fuchsian singularity of the OUT-system (16), with residue —A^/X. Therefore, we can 
always find a fundamental matrix solution with behavior: 



tf.Match 

^OUT ~ 



I + O 



OO. 



This solution matches with Also A = 1 is a fuchsian singularity of (16). Therefore, we have: 

^atch =jOUT^[ I + ( X _ 1 ^( X _ 1 ft* a ( X _ 1 )Ri C OUT i A 1; 

Here C° UT is a suitable connection matrix. ipi UT (x) is the matrix that diagonalizes the leading terms 
of Ai(x). Therefore, -01 {x) ~ ^P\ UT i x ) f° r x —¥ 0. As a consequence of isomonodromicity, i?i is the 
same of 

As a consequence of the matching <3> o ^QUT h ^ the monodromy of ^ at A = 1 is: 



OUT 



Mi = M° VT = CT 1 exp{i7r0iCT 3 } exp{27rii2i}Ci, with d = C{ 

We finally need an invertible connection matrix Cout to connect ^Q^ h with the solution ^out 
which was matching with ^/jv- Namely 

W Ot/T — ^OUT^OUT- 



(*) Matching * «-> $ 



7 AT- 



As a consequence of the matching <h> ^qut > we have to choose the IN-solution of (17) which 



matches with *^ c \ This is 



,t, Match . ,T f >r 



Now, A = 0,x are fuchsian singularities of the IN-system. Therefore: 

^ JV (x)[/ + 0(A)] A^A^Cft* A^O; 
^ Ar (x)[/ + 0(A-x)] (A-x)^ ff3 (A-x)«*Ci w , A -> x; 



vr.Match 



The above hold for fixed small i^O. Here Cq^ and are suitable connection matrices. i/jq N (x) and 
ipx{x) IN are diagonalizing matrices of the leading terms of Aq{x) and A x {x). For x — > they match with 
if>o(x) and ^ x (x) of ^ in (55). On the other hand, as a consequence of isomonodromicity, the matrices 
i?o and R x are the same of The above ^f* tch has the same behavior of "J at A — > and A — > x; 
moreover, it is an approximation of * for x small. The matrices Cq N , are independent of x. So, 
the matching o ^/jv is realized and the connection matrices Co and C x coincide with Cq N , 
respectively. As a result, we obtain the monodromy matrices for *: 



M = M /Ar = GT 1 exp{iTT0 O (T 3 } exp{27rii?o}C , C = C, 



o > 



M x = M K /Ar - C^ 1 exp{i7r0 x <7 3 } exp{27rifl x }C x , C x = C™. 

The matrices Ci UT ', Cout, Cq N and C x w can be computed because we are able to solve exactly (i.e. 
in terms of special functions) the reduced systems OUT and IN. Details are technically complicated, and 
the reader is refcrcd to [10], [11]. 
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10 Conclusions 



We have given the table of critical behaviors at x — of PVI-transcendents, obtained with the method of 
monodromy preserving deformations (essentially, the Jimbo's procedure [17] and the matching procedure 
[10]). From the table, all the behaviors at x = 1 and oo can be computed via a birational transformation 
(12) and (13). We have explained in details the meaning of the behaviors, and reviewed how they were 
obtained. We have also given the parameterization in terms of monodromy data of the basic critical 
behaviors. The parameterization allows to solve the connection problem. We have explained that the 
table of critical behaviors is essentially complete. We have shown that the expansions obtained by power 
geometry in [6] coincide with the behaviors obtained by the monodromy preserving deformation method. 

This paper is aimed at providing a synthetic (and hopefully useful) tabulation of both critical be- 
haviors and parametrization formulas, which so far has been missing. 

Acknowledgments: I thank the Korea Institute of Advanced Study (KIAS, Seoul 130-722, Republic 
of Korea) for providing computing resources (Abacus System) for this work. 
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